Exotic Superconducting Properties in the Electron-Hole Compensated Heavy Fermion 

'SemimetaP URu 2 Si 2 
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We show that the charge and thermal transport measurements on ultraclean crystals of URu2Si2 
reveal a number of unprecedented superconducting properties. The uniqueness is best highlighted 
by the peculiar field dependence of thermal conductivity including the first order transition at 
H C 2 with a reduction of entropy flow. This is a consequence of multi-band superconductivity with 
compensated electronic structure in the hidden order state of this system. We provide strong 
evidence for a new type of unconventional superconductivity with two distinct gaps having different 
nodal topology. 



The heavy- Fermion compound URu2Si2 has mys- 
tified researchers since the superconducting state 
(T c = 1.5 K) is embedded within the "hidden order" 
phase (T h = 17.5 K) 0, H 0|. Although several ex- 
otic order parameters have been proposed for the hid- 
den order phase [1], it is not identified yet. According 
to several experimental observations, most of the carri- 
ers disappear below T/j resulting in a density one order 
of magnitude smaller than in other heavy- Fermion super- 
conductors @, @j 0] • Superconductivity with such a low 
density is rcmarkablcsincc the supcrfluid density is very 
low in some way reminiscent of underdoped cuprates. 
Moreover, pressure studies reveal that the superconduc- 
tivity coexists with the hidden order but is suppressed 
by antiferromagnetic ordering Q. 

In this Letter, using ultraclean single crystals, we re- 
port various anomalous superconducting properties in 
URu2Si2- We show that a peculiar electronic structure 
appearing below the hidden order transition provides an 
intriguing stage on which a new type of unconventional 
superconducting state appears. 

Single crystals of URu2Si2 were grown by the Czochral- 
ski pulling method in a tetra-arc furnace. The well- 
defined superconducting transition was confirmed by the 
specific heat measurements. The thermal conductivity 
k was measured using a standard four-wire steady state 
method along the a-axis (heat current q\\ a). The contact 
resistance at low temperatures is less than 10 mfl. 

We first discuss the electronic structure below Th ■ Fig- 
ure 1 shows the temperature dependence of the resistivity 
p along the a-axis and Hall coefficient Rh (solid circles) 
defined as Rh = -§W- at H — > T for H\\ c in the tetrag- 
onal crystal structure. In zero field, p depends on T as 
p = po + AT 2 below 6 K down to T c . The exceptionally 
low residual resistivity po = 0.48 fj,Q cm and large resid- 
ual resistivity ratio RRR = 670 attest the highest crystal 
quality currently achievable. Rh exhibits an eight-fold 



increase below Th and attains a T-indcpcndcnt value at 
low temperatures, associated with a strong reduction of 
the carrier density. Most remarkably, the magnctorcsis- 
tance (MR) increases with decreasing temperature and 
becomes extremely large at the lowest temperatures. The 
inset of Fig. 1 shows the MR Ap/p = (p(H) — po)/po 
extrapolated for T — > K as a function of H 2 . The MR 
exhibits a nearly perfect independence without satura- 
tion up to 10 T. At poH = 10 T, Ap/po reaches as large 
as 300. We note that the observed MR is much larger 
than in any known superconducting material. 

The combination of the peculiar electronic structure 
of URu2Si2, i.e. high purity, the low carrier density and 
essentially equal number of electrons and holes (compen- 
sation), n e = Uh, gives rise to this exceptionally large 
MR [9(. The observed quadratic field dependence of the 
MR without saturating behavior demonstrates the nearly 
perfect compensation, is given by 



Ap/po = (uj e c T e ) ■ (uj'cT h ) 



and the Hall resistivity as 
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where uj c = eH/m is the cyclotron frequency for carriers 
of the mass m and the scattering time r Q . The suffixes 
e and h denote 'electron' and 'hole', respectively. From 
Eq. (1), we estimate {oj e c r e ) ■ (u%T h ) ~ 25 at iJ c c 2 = 2.8 T, 
where H^ 2 is the upper critical field for H\\ c. We note 
that this large w c r-valuc sets a record of type-II super- 
conductors. The best-known systems exhibiting a similar 
behavior are very pure semimetals, such as bismuth and 
graphite [l(j. The T-independent positive Rh at low 
temperatures together with nearly iJ-linear p xy up to 
10 T (not shown) indicates that wJt/, ^> uj^.t c , or r e and 
Th have the same T-dependence. In cither case, 1/Rne 
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FIG. 1: (color online). Main panel: Temperature dependence 
of p (solid lines) along the a-axis in zero and finite magnetic 
fields (H\\ c) and Hall coefficient Rh (solid circles). Inset: 
The MR at T -> K as a function of H 2 . 



gives an upper limit of the hole number. The number 
of holes estimated from Rh = 2.3 x 10~ 8 m 3 /C at low 
temperatures, where the influence of the skew scattering 
can be neglected, is nt % 2.6 x 10 26 m -3 (corresponding 
to 0.021 holes per U-atom). Such a low carrier density 
implies Fermi surfaces occupying only a small part of 
the Brillouin zone. According to de Haas- van Alphcn 
(dHvA) experiments [ll|, [T3] , one hole band with nearly 
spherical shape is reported. This band with Fermi mo- 
mentum kp ~ 1.8 x 10 9 m _1 contains ~ 0.017 holes/U. 
This number gives a lower limit of hole number, because 
of the presence of very small hole pockets. Thus we ob- 
tain 0.017 < n h < 0.021 holcs/U. 

The compensation necessarily indicates the presence 
of electron band, although it is not confirmed by the 
dHvA experiments lll|, [l2[ . The magnitude of the elec- 
tronic specific heat [2[ is five times larger than the esti- 
mation taking only the hole band assuming = 13mo 
This indicates that the electron mass is much 
heavier than hole mass. Moreover, the observed large 



anisotropies of the MR and the upper critical fields [11 1 
indicates that the shape of the heavy electron band is el- 
liptical [r| • We note that heavy electrons are consistent 
with the Hall effect which is dominated by light holes. 

We next address the superconducting state. The inset 
of Fig. 2 shows the T-dcpcndence of the thermal con- 
ductivity divided by T, k/T. The Wiedemann- Franz 
(WF) ratio L = (n/T)p at T c is (1.00 ± 0.02)i , where 
L Q = 2.44 x 10~ 8 flW/K is the Sommerfeld value, in- 
dicating that the electron contribution dominates over 
the phonon contribution. Below T c , k/T displays a steep 
increase and exhibits a pronounced maximum at 0.5 K, 
which arises from the strong suppression of the quasipar- 
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FIG. 2: (color online). Main panel: k/T plotted as a function 
of T 2 at very low temperatures. The data in magnetic fields 
(if || c) are obtained under field cooling conditions. The lines 
are guides for eyes. Inset: k/T vs T at high temperatures. 



ticle (QP) scattering rate due to the formation of a gap, 
as observed in several strongly correlated systems [14| . 

Figure 2 depicts k/T plotted as a function of T 2 at 
very low temperature. In zero field, k/T is well fitted 
by k/T — kqq/T + bT 2 , where b is a constant. There 
is a fundamental difference in the low-temperature heat 
transport between superconductors with isotropic and 
anisotropic gaps. In the latter, the QP heat conduction 
is entirely governed by nodal topology (lo| . As shown in 
Fig. 2, the presence of a residual value in k/T at T — > K, 
Koo/T, is clearly resolved in zero field. This indicates the 
presence of the residual normal fluid, which can be at- 
tributed to the presence of nodes in the gap function. A 
more detailed analysis will be given later. 

The iJ-dcpcndence of k(H)/T obtained by the extrap- 
olation to T — > K at each field value (see Fig. 2) is quite 
unusual, as shown in Figs. 3 (a) and (b). There are three 
characteristic fields denoted as H s ~ 0.4 T, H^ 2 ~ 2.8 T 
and H^ 2 — 12 T, where ff™ 2 are the upper critical field 
for H\\ a. For both H\\ a and Jf|| c, k{H)/T exhibits 
a pronounced increase below ~ H s (Fig. 3(b)), followed 
by a plateau- like behavior above ~ H s . Interestingly, in 
spite of the large electronic anisotropy, k(H)/T is nearly 
isotropic with respect to the field direction at low fields, 
as shown in Fig. 3(b). For H\\ c, k(H)/T exhibits a jump 
to a very low value, which is nearly 1/4 of k(H)/T just 
below 11^2, in a narrow field window (AH < 0.2 T) at 
H^2- The dashed blue and dotted red lines represent the 
expected k(H)/T by the WF-law for H\\ c and a, 
respectively, assuming the quadratic dependent MR (see 
the inset of Fig. 1). Above H^ 2 , k(H)/T decreases rapidly 
with H related to the remarkable MR, following the WF 
law k(H)/T = L Q /p(H). On the other hand, for H || a, 
k(H)/T increases steeply with H with downward curva- 
ture over a wide range of fields approaching H® 2 - Near 
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FIG. 3: (color online), (a) The field dependence of k/T for 
Jf|| c (blue circles) and H || a (red squares) as a function of 
H/Hc 2 and H / H^ 2 , respectively. The dashed blue and dotted 
red lines represent the expected k(H)/T by the WF-law for 
(H\\ c) and (H\\ a), respectively, (b) n(H)/T at low fields 
as a function of \fH. The characteristic field, H s =0.4 T is 
indicated by the arrow, (c) The field derivative of k(H)/T 
for Jf || c plotted as a function of H. (d) The H — T phase 
diagram for H || c. The transition between superconducting 
(SC) and normal (N) states is of first order below ~ 0.4 K. 



H®2, k(H)/T turns to decrease. It is natural to assume 
that above our experimental field limit (12.4 T), k(H)/T 
drops to a lower value obtained by the WF law, similar 
to the case for if|| c. 

To our knowledge, the remarkable step-like reduction 
of k(H)/T at upper critical fields is unprecedented in 
type-II superconductors. Since the thermal conductivity 
is related to an entropy flow, its jump (irrespective of its 
directions) immediately indicates the occurrence of the 
first order phase transition, as reported in CeCoIns [lij ]. 
To determine the end point of the first order transition, 
we plot d(-^)/dH vs. H at temperatures below 0.5 K for 
H\\ c in Fig. 3(c). A distinct dip structure can be seen 
below 0.4 K but is absent at 0.5 K, indicating that the 
end point is at ~ 0.5 K. Figure 3(d) illustrates the H — T 
phase diagram for H\\ c. Since the Maki parameter of 
URu2Si2 is more than unity for both field directions, the 
origin of the first-order transition is most likely due to the 
Pauli paramagnetic limiting. This suggests that URii2Si2 
realizes an even-parity spin-singlet pairing state. 

Usually, the heat conduction just below the upper crit- 
ical field H C 2 is smaller than that in the normal state, 



owing to the reduction of number of QPs carrying heat 
and the enhancement of the scattering rate by vortices. 
Therefore the observed field induced transition at H C 2 
of URu2Si2 with an enhanced entropy flow in the super- 
conducting state compared to the normal state is highly 
unusual. We point out that the compensation effect pro- 
vides an explanation. In the normal state, the compen- 
sation leads to an extremely large MR and hence large 
reduction of the thermal conductivity with increasing H . 
However, in the superconducting mixed phase, the com- 
pensation of hole- and electron-like QPs is violated be- 
cause the size and the nodal structure of the gaps on the 
two Fermi surfaces are very different. This implies that 
the heat transport is drastically changed at H C 2- 

The initial steep increase and subsequent plateau be- 
havior of the thermal conductivity shown in Figs. 3(a) 
and (b) indicate that a substantial portion of QPs is al- 
ready restored at H s much below the upper critical fields. 
Such a field dependence bears a strong resemblance to 
those of MgB 2 [13 and PrOs4Sbi2 [18| . providing evi- 
dence for multi-band superconductivity in URii2Si2; H s 
is interpreted as a virtual upper critical field that con- 
trols the field dependence of the smaller gap. The multi- 
band superconductivity is, of course, compatible with the 
compensation effect originating from the presence of two 
bands. It is natural to assume that the smaller gap origi- 
nates from the light hole band. In addition, the spherical 
structure of hole band is consistent with nearly isotropic 
k(H)/T below H s . 

The steep increase of n(H)/T at low fields is in strong 
contrast to the behavior of fully gapped superconductors, 
where all the QP states are bound to vortex cores and, 
therefore, the applied magnetic field hardly affects the 
thermal conduction except in the vicinity of H C 2- By con- 
trast, in superconductors with nodes, the heat transport 
is dominated by the delocalizcd QPs and the Dopplcr 
shift of the QP spectrum induces a remarkable field de- 
pendence of k(H)/T. For line nodes, the QP density 
of states N(0) increases with \/H [l5j. As shown in 
Fig. 3(b), k(H)/T cx N(0) for both field directions is 
nearly proportional to \Jli up to ~ H s . The presence of 
line nodes is further supported by a finite k o /T term in 
zero field. For line nodes, kqo/T is given by ^ ~ ^ • ^ , 
where £ a b and £ e are the coherence length and mean free 
path within the a6-plane, respectively [l9| . To estimate 
£ Q i,, we use H s instead of H^ 2 , assuming the line nodes in 
the smaller gap. Using £ e ~ 1.2 (jxn in our crystal [20l ]. 
koo/T is estimated to be 0.15 W/K 2 m. This value is in 
the same order with the observed kqo/T = 0.38 W/K 2 m. 
These results suggest the presence of line nodes in the 
gap of the light hole band. 

We next discuss the gap structure of heavy electron 
band. This can be inferred from the thermal conductiv- 
ity well above H s where essentially all QPs of the light 
band (small gap) contribute to the heat transport. The 
plateau-like behavior of k(H)/T up to H£ 2 for H\\ c in- 
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FIG. 4: (color online). The schematic figure of the Fermi 
surface (opaque) and superconducting gap structure (trans- 
parent) inferred from the present experiments [2lj . 

dicates that the QPs of the heavier band do not expe- 
rience the Dopplcr shift. On the other hand, for H\\ a, 
the steep increase of k(H)/T with downward curvature 
above ~ 0.2iJ" 2 up to H " 2 can be interpreted as a Doppler 
shift effect for the QPs due to nodes. Note that in fully 
gapped superconductors, k(H)/T increases steeply only 
just below the upper critical field with an upward cur- 
vature. For complete understanding of the multi-band 
effects, the quantitative description for the crossover be- 
havior from plateau to Doppler regimes for H\\ a is re- 
quired. However, the fact that the Doppler shift occurs 
only in parallel field suggests that it originates from point 
nodes in the gap along the c-axis as these would not yield 
a Dopplcr shift for fields along the c-axis. Line nodes in 
contrast would give rise to a Dopplcr shift for any field 
direction (IBj . In this way we distinguish the node topol- 
ogy in the gaps on the two bands. 

The present experiments strongly suggest (1) line node 
in the spherical light hole band, (2) point nodes in the 
elliptical heavy electron band (see Fig. 3(a)), and (3) 
spin-singlet pairing. Although the detailed Fermi sur- 
face topology has not been fully determined, it is tempt- 
ing to speculate on the possible gap symmetry which 
is consistent with the experimental observations. The 
classification scheme of even parity Cooper pairs allows 
for one type of state which has simultaneously (sym- 
metry protected) line and point nodes. This is a two- 
component order parameter combined in a time reversal 
symmetry violating combination of the basic <i-wave form 
Afc oc k z (k x ± iky). As the nodal structure is different 
on the two types of Fermi surfaces, it is necessary to ac- 
count for the behavior of gap in the entire Brillouin zone. 
For the body-centered tetragonal crystal lattice the basic 
form of the gap function with the correct periodicity in 
fc-space is given by 

a a fez ( fex fey fex fey \ / \ 

A = An sin — c sin -a + i sin -a , (3) 

2 V 2 2 J w 

resulting from (nearest neighbor) intcrlayer pairing. A 



schematic form of the gap structure and the Fermi pock- 
ets in the Brillouin zone is depicted in Fig. 4 [21]. The 
line nodes lie on the zone boundary on the hole-like light 
band, while both Fermi surfaces have point nodes along 
the c-axis. In avoiding a Fermi surface centered at the 
T-point the line node for k z = does not play a role. 
This superconducting phase possesses a two-component 
order parameter which may give rise to interesting fur- 
ther features beyond the nodal properties of the QP gap. 

In summary, the charge and thermal transport mea- 
surements on ultraclean URu2Si2 reveal a number of un- 
precedented superconducting properties. We provide ev- 
idence for superconductivity with two distinct gaps hav- 
ing different nodal topology. This intriguing supercon- 
ducting state adds a unique example to the list of uncon- 
ventional superconductors. 
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specific heat study which also suggests the point nodes 
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